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Abstract. Chcmotaxis is typically modeled in the context of cellular motion towards 
a static, exogenous source of chemoattractant. Here, we propose a time-dependent 
mechanism of chemotaxis in which a self-propelled particle {e.g., a cell) releases a 
chemical that diffuses to fixed particles (targets) and signals the production of a second 
chemical by these targets. The particle then moves up concentration gradients of this 
second chemical, analogous to diffusive echolocation. When one target is present, we 
describe probe release strategies that optimize travel of the cell to the target. In 
the presence of multiple targets, the one selected by the cell depends on the strength 
and, interestingly, on the frequency of probe chemical release. Although involving an 
additional chemical signaling step, our chemical "pinging" hypothesis allows for greater 
flexibility in regulating target selection, as seen in a number of physical or biological 
realizations. 
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1. Introduction 

Organisms that employ chemical signaling for functions such as antimicrobial defense 
mechanisms [1, 2, 3] and nutrient uptake [4, 5] must coordinate the influences of a 
complicated spatio-temporal mixture of signaling molecules emanating from possibly 
many sources. A fundamental problem in chemotaxis is how the cell determines a 
strategy to best select a specific target from many? This problem arises for many 
chemotactic organisms, such as E. coli [6], Myxococcus xanthus [7], and Dictyostelium 
discoideum [8], and poses a challenging theoretical question. 

Here, we propose a biologically plausible mechanism involving active sensing, in 
which chemical signaling is initiated by a chemical prompt from the chemotactic cell. 
Potential biological manifestations of active sensing include "diffusion sensing" [9] and 
cancer cell chemotaxis [10]. Diffusion sensing is a particular variant of quorum sensing 
that involves release of a metabolically inexpensive compound, inducing nearby cells to 
emit the main signal back to the first cell indicating to it that that other chemically 
responsive targets are in the vicinity [9]. Cancer cells also exploit this type of indirect 
sensing by producing excess amounts of matrix metalloproteinases (MMP) which cleave 
substrates (such as growth factors and laminins) bound to the extracellular matrix 
(ECM) [10]. These cleavage products diffuse back to the cancer cells, inducing their 
chemotaxis. Dynamic multistep chemotaxis mechanisms may also arise in the mating of 
yeast, where each sex of yeast emits its own pheromone that up-regulates gene expression 
of the complementary pheromone in the opposite sex [11]. 

In our model of active sensing, moving chemical sources and targets interact 
through the time-dependent diffusion of signahng molecules Therefore, the timing of 
the release and detection of chemoattractants will also be important in the overall 
chemotactic process [12]. Indeed, there is evidence that the probe emission of MMP's 
is transcriptionally regulated and is time-dependent [13]. Experiments have also 
demonstrated the existence of robust and tunable oscillations in transcription in E. coli 
[14] and mammalian cells [15], and cAMP release in Dictyostelium discoideum [8, 16]. 
Another important virtue of the proposed dynamic sensing mechanism is that it allows 
cells to detect local, transient chemoattractant gradients in complex media where steady- 
state gradients cannot be sustained. For example, branched, "dead-end" volumes with 
impenetrable boundaries cannot support a steady-state chemical gradient, but can allow 
a transient gradient. Therefore, chemotaxis under time-dependent, and in particular, 
time periodic conditions are realizable systems for further exploration. 

Basic models for these potentially novel time-periodic chemotactic systems are 
currently lacking. Here, we distill the extracellular components of a dynamic multistep 
chemotaxis mechanism into an essential physical model that describes gradient-guided 
cellular motion towards chemically responsive targets. Classic chemotaxis models 
by Keller and Segel [17], and their extensions [12, 18] consider passive sources of a 
chemoattractant to which a cell responds, or consumption of a passive chemical [19]. 
Contrary to the classic model of chemotaxis where the cell moves along a static chemical 
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Figure 1. A schematic of a motile cell releasing probe chemical a (blue) which gets 
converted at fixed targets to a chemoattractant b (red) which the cell detects and uses 
to move. For the two target case shown here, the initially farther Target 2, has a larger 
a — > 6 conversion rate (^2 > fci)- 

gradient of nutrients already present in the environment, we consider the following 
scenario. Initially, a cell sends out a chemical probe signal (species a), which diffuses 
to stationary targets. Upon contact with the targets, the probe chemical a induces 
the target to release a different chemical (species b), which diffuses back to the cell. 
In response, the cell moves up the chemical gradient of the chemoattractant b. A 
schematic of our proposed "chemolocation" mechanism, a diffusive analogue of sonar, 
is shown in Fig. 1, where one cell and two targets are depicted. Here, we assume that 
chemoattractant b is different from probe chemical a ("paracrine" signaling) and that 
the cell uses only b to guide its motion towards the target (s). 

2. Mathematical Model 

Denote the concentrations at spatial position r of probe chemical a emitted by the cell, 
and chemoattractant b produced by the fixed targets by na(r, t) and nb{r, t), respectively. 
The production of chemoattractant b by the targets may be initiated by binding of probe 
a to receptors on the targets. Although the cell needs to detect spatial gradients in rib, 
we assume for simplicity that both the cell and the target can be treated as point 
particles when considering the diffusive dynamics of Ua and Ub- 
The governing equations in our model are 

hair, t) = DaV^Ua - fiaUa + F{t)6{r - R(t)), (l) 

hb{r,t) = DbV^Ub - fibUb + ^5(r - Rj)Kj[na{r,t - t'),t] , (2) 

j 
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m= f 7(t')Vt/K(r,t-t')]dtl=R(,) (3) 

where 5{-) is the Dirac delta function, R(t) is the position of the moving cell, Rj are the 
fixed target positions. Da, are the uniform probe and chemoattractant diffusivities, 
and fia, fib their uniform degradation rates. In Eq. 1, F{t) represents the time-dependent 
emission of probe chemical a by the cell at position R(t), while in Eq. 2 the term 
^(r ~ [^a(r,t — ^0)^] represents the total production of b from a by all the 

targets j, each at fixed position Rj. Note that for systems in dimension o? < 2, we 
must assume Hail^h > if F{t — > oo) > in order for the chemical concentrations to 
remain bounded. The functional Kj [na(Rj, t — t'),t] embodies all chemical steps in the 
production of chemoattractant b in response to contact with probe chemical a. The 
production of chemoattractant b occurs with a delay t' , (the time taken for a single or 
multi-step reaction along the signaling pathway) after the targets are exposed to the 
probe a at a time t — t' [4, 20]. We will henceforth assume the simplest model for the 
kernel Kj [naCR-jji — t'),t] = kjna(R.j,t), representing instantaneous production, with 
rate kj, of chemoattractant proportional to the concentration of probe. 

Equation 3 describes the motion of the cell in an effective time-dependent potential 
U [nb{'R{t),t)] generated by the dynamics of the chemoattractant nb{'R{t),t). The 
functional U [nb{r, t)] could be a nonlinear function of ^^(r, t), such as one with threshold 
and saturation like U oc n"(r, t) / (const. -|- n"(r, t)), where a is the Hill coefficient. This 
form may be more appropriate when cooperative binding of many chemoattractant 
molecules is required to trigger cellular migration. The response 7 itself may also be 
time-dependent. Also, note that not only might the production of chemoattractant b 
encounter a delay, the cell response embodied by 7 might be also delayed. Although such 
delayed and nonlinear responses may result in intrinsically rich signaling behavior, for 
the sake of simplicity, and in order to analyse our results with as few free parameters as 
possible, we assume the "force" on the cell is proportional to the local chemoattractant 
gradient i.e., 

R(t)=7Vn6(r,t)|r=R(,), (4) 

Stochastic effects due to low Ua, Ub concentrations or other random effects within the 
signaling process [21] can also be easily incorporated component by either considering 
a randomly varying j(t) and/or by adding a random noise to Eq. 3 or 4. The additive 
noise is equivalent to endowing the cell with Brownian motion. Nonetheless, we show 
that even in the ideal case of perfect signaling and constant 7, novel target selection 
phenomena arise. 

In the following analysis, we define dimensionless parameters by measuring length 
in units of the initial separation between the probe and its farthest target, and time 
in units of {R^Y/Da- In spatial dimension d, the dimensionless equations are identical 
in form to Eqs. 1-3 except with Da = 1, and the terms Db, fia,b, kj, 7, and F(t) replaced 
by Db/Da, fia,bRl/Da, kj/{DaRt^), ^/{DaRf), and RlF{t)/Da, respectively. 
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The solutions to Eqs. 1 and 2 can be solved in terms of Green's functions and the 
cell position R(t). Upon substitution of % into Eq. 4, we find a self-consistent nonlinear 
equation for the cell position 

R(t) =-f^kj [ dt' I dt"F{t")Ga{R{t") - RjS t' - t")Hfe(R(t) - Rj; t - t'), (5) 
j Jo Jo 

where Ga{r;t) = (47rt)-'^/2g-r2/4tg-/.„t^ and 

ll,{R{t) - Rj, t-t') = VG,{v - R,-; t - t')l^^^^~^ 

_ 27r(R,- - R(t))e-|i^W-i^^lV[4i?(*-*')]e-M6{t-t') (6) 
~ [AnDit - t')]^/2+i ■ 

Although the model equations are linear in Ua, rif,, the moving source of probe chemical 
renders the problem intrinsically nonlinear and not amenable to analytic treatment. 
Since bounds and analytic expressions for n},{r, t) and R{t) can be found only in special 
cases (such as F{t) oc S{t)), we will solve for the cell trajectory by either numerically 
integrating Eq. 5, or by directly numerically solving Eqs. 1-3 on a fixed lattice using 
a stable backward-time, central space scheme with step sizes Ax = At = 10"'^. in this 
case, the system boundary is chosen to be far enough away from the targets as to be 
irrelevant. We have verified that our results do not depend on the numerical approach 
employed. 

3. Results and Discussion 

We shall study our model predominately by solving either Eqs. 1-3 or Eq. 5 numerically, 
and exploring the qualitative features of the chemolocation mechanism. However, 
in certain physical limits, we find analytical relationships useful in describing target 
selection. 

3.1. Single Target 

In biological media or in laboratory realizations [22], diffusion often occurs in confined 
or ramified geometries. For chemotaxis across capillaries, or across percolating paths, 
the effective dimensionality d of the diffusion process may be smaller than the spatial 
dimension. For simplicity, we first explore the qualitative behavior of the one- 
dimensional {d — 1 and R{t) = X{t)) version of our model with initial condition 
na{x, 0) = ni,{x, 0) = -'^(O) = 0, and 7 = 1. As a demonstration of the chemolocation 
mechanism, consider a cell moving towards a single target under different probe release 
protocols F{t). Different strategies of probe release qualitatively infiuence the ability of 
the cell to reach the target. 

Figure 2 shows the trajectories and maximum distance traveled by the cell when 
a single (5-function probe is emitted at t = 0. Here, and in the rest of the paper, the 
"5-function" is approximated by a narrow square pulse release of duration dt = 0.1 
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Figure 2. Chemolocation to a single target using a single 5-function probe release. 

(a) Motion of the probe towards a single target in d = 1 using three different 
(5-function release intensities. Delta-functions were approximated with releases of 
duration dt = 0.1. Parameters used were D — j = k — 1 and /ia = ^J■b = 0.1. 

(b) Critical signaling strength ^kQ as a function of chemical decay rates fi = ^ fJ-b- 
Results from numerical solutions and the upper bound Eq. 8 are shown for both d = 1 
and d — 2. 



and intensity Q. The cell starts to move only after a short delay during which some 
of the probe has reached the target, and the converted chemoattractant has diffused 
back to the cell. For a single impulse release of probe a, the velocity of the cell towards 
the target is initially high but eventually goes to zero since the system runs out of 
the chemoattractant b once the single pulse of probe a has dissipated. Thus, for a 
modest {e.g. Q = 1,2) single-pulse release in Fig. 2(a), the cell moves only part of 
the distance to the target. In the low mobility (small 7) limit, an approximation to 
the total travel distance X^o after probe release in a single pulse can be obtained by 
setting R(t") = R'(^) = in Eq. 5. In one-dimension {d = 1), integrating the left-and- 
side of Eq. 5 yields X(t)dt = X^q. Explicit evaluation of the time integral of the 
right-hand-side of Eq. 5 shows that 

V r J I Jl 

When there is only a single target at Xi = 1, Eq. 7 gives a lower bound for X^o, 
implying that 

^kQ > AD^e^^y^ {d = 1) (8) 

is a sufficient condition for the cell to reach the target. The analogous sufficient 
conditions for insuring the cell reaches a single target in c? = 2 and d = 3, are 

ikQ > (rf = 2), (9) 
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and 



7^<3 > 



((i = 3), 



(10) 



respectively. These "critical" values of the signaUing strength ^kQ represent the 
analytically obtainable lowest values above which the target is acquired with certainty. 
For d — 2, 3, these conditions derived from the R(t") = R'(^) = approximation are 
independent of the relative diffusivity D. In higher dimensions, the spreading of probe 
and chemoattractant concentrations in more severe so that under the approximation 
R(t") = R(t) ~ 0, the stationary cell experiences a significantly diminished gradient 
compared to that which it experiences if it had moved. Chemical decay /j, also 
diminishes the signal. Therefore, the signalling strength required for reaching the target 
is significantly increased in this "adiabatic" approximation, especially for large d or /i. 
Fig. 2(b) plots the critical values (for d—1 and d — 2, d — 3 not shown) of the signaling 
strength jkQ, from both simulations and Eq. 8, as functions of = = /ib- Indeed, 
the sufficient condition Eq. 8 for = 1 provides the tightest bound. 

The likelihood of arrival to a single target can be enhanced not only by increasing 
the signaling strength '~fkQ, but also by releasing multiple pulses (as will be shown in 
Fig. 4(b)) and by releasing probe more slowly. Consider the case where the cell contains 
a fixed amount of probe chemical a. How should the cell release this fixed amount of 
probe to best reach the target? Suppose the release occurs in a single pulse of duration r 
and amplitude Fq — F{t < r) such that the total amount Q — Fqt of chemical released 
is constant. Fig. 3(a) shows trajectories for a cell that releases a fixed amount, Q, of 
chemical a at a constant rate Fq over different lengths of time r = Q/Fq. For Fq = 20 
and r = 1, the large magnitude, short duration release allows the cell to travel only 
approximately 90% of the way to the target. In such cases, the cell's velocity may reach 
a high value; however, the cell docs not reach the target before all chemical signals have 
dissipated. For lower intensity, but longer duration releases, the cell is better able to 
reach the target. Quahtatively, this can be understood by noting that a single 5-function 
release gives a lower bound on the distance traveled. If the probe is released more 
slowly, the cell moves shghtly closer to the target, amphfying the effect of the probe 
chemical that is released at later times, because this portion can reach the cell more 
quickly reducing its decay. 

In the limit r — )■ oo, with fixed Q = FqT, we can find an upper bound X™^^ 
for the distance travelled to a single target by integrating Eq. 7 and summing the 
total distances travelled for each independent increment dQ = -fodr of probe release. 
Assuming = /x^ = A'-, -D = 1, and Xi = 1 for simplicity, we find for ci = 1, 



valid for < '~fkQ < 2D[e'^^ — 1). Therefore, even for infinitely slow release of a fixed 
amount Q of probe chemical, a necessary condition -/kQ > 2D{e^^ - 1) remams, even 
if the cell can detect arbitrarily low concentration gradients of chemoattractant. 




(11) 
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Even if the cell has sufficient probe to reach the target, if the probe is released very 
slowly, the time required for the cell to reach the target may be large. Fig. 3(b) plots 
the time t* it takes for the cell to reach the target, as a function of the pulse duration 
r. For large signaling strength 'ykQ, the cell reaches its target for all r, and reaches 
its target most quickly for very short pulse durations, r. When the signaling strength 
jkQ is decreased, the cell reaches its target only when r is greater than some critical 
value, denoted by the vertical asymptotes (dashed lines) in Fig. 3(b). Additionally, 
there is a value of r which minimizes the arrival time t*, and approximately satisfies the 
condition t ^ t*, indicating that the cell reaches its target in the shortest time when 
probe chemical a is released over a period r nearly equal to its travel time. If we think 
of probe a as an effective "fuel" that drives the motion of the cell, it is not ideal to have 
r <C t*, because if chemical a stops being released well before the cell reaches its target, 
the cell will have to rely on a residual, decaying concentration field b to reach its target. 
On the other hand, if r 3> t*, the cell will continue to release chemical a after reaching 
its target, when a faster release of a would have allowed the cell to reach its target more 
quickly. 
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Figure 3. Chemolocation to a single target using a probe pulse of duration r. (a) 
Trajectories of a cell to a single target when a fixed, total amount of probe Q is released 
at a constant rate, Fq for a finite time, r. (b) Arrival times t* of a cell to a single 
target are plotted as a function of pulse duration r. Within each curve, the total 
amount of probe released, Q was held constant. In both plots, parameters used were 
L) = 7 = fc = = Mb = 1- 



3.2. Multiple Targets 

We now illustrate the mechanism of frequency and response strength-dependent target 
selection among multiple targets. First consider the one- dimensional case with a cell and 
two targets. The closer target (at Xi = —0.25) is assigned a probe-to-chemoattractant 
production rate ki = 0.5, while the farther target at X2 = 1 has a larger production 
rate /c2 = 1- Fig. 4(a) shows that when the release of probe chemical is in the form of a 
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Heaviside function, target selection can be controlled by the cell's response, 7. For the 
parameters chosen, if 7 = 0.1, the far target is chosen. As the strength of the response 
is increased beyond approximately 0.225, the near target is selected by our chemotactic 
mechanism. When the cell's response, 7 is large, the cell is initially pulled strongly 
towards the near target. Because the distance to the far target increases substantially 
before the signal from this target reaches the cell, the signal strength diminishes and 
is insufficient to pull the cell back towards the far target. When the cell's response is 
small, the cell is unable to move much towards the near target before the stronger signal 
from the far target reaches the cell, and the cell ultimately gets pulled to the far target. 



9(r) probe release pulsed probe release target selection 
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Figure 4. Target selection as a function of chemotactic response 7 and probe release 
frequency 1/T. (a) Response-dependent target selection showing a transition from 
selecting the further, stronger, target to selecting the weaker, nearer target for a 
Heaviside release protocol F{t) — 9{t) and a varying response 7. When 7 > 0.235, the 
near target is selected, while the far target is selected for smaller 7. (b) Frequency- 
dependent target selection, with release protocol F{t) — T^-S{t — iT) and varied 
response strength 7. The weaker nearer target is selected for constant chemical release 
F{t) = 9{t) and very high frequency pulsed release F{t) = 0.2^. (5(i — 0.2i), while 
the further stronger target is selected for pulsed release at intermediate frequencies 
(corresponding to interpulse intervals T — 3,7). These trajectories may first arrive 
at the nearer target depending upon probe chemical release amplitude, (c) Under the 
release function F{t) = Tj^i ^{t ~ ^T), regions of 7 — T space where the cell reaches 
the nearer (and weaker), or farther (and stronger) target are shown. When T > A, 
the cell always reaches the far target. For smaller values of T, the cell reaches the 
far target for small and large values of 7, and reaches the near target for intermediate 
values of 7. For all plots, the target strengths are ki = 0.5, k2 — I, the decay rates are 
Ma — — 0.1 and D = 1. The initial separations between the cell and the targets are 
\Xi\ = 0.25 and \X2\ = 1 respectively. 

Target selection can depend not only on probe release intensity and response 
strength, but is also modulated release frequency. Fig. 4(b) illustrates target selection 
when the probe chemical is released either as a Heaviside function 6{t), or as a series 
of pulses F{t) = T'Y^^5{t — iT) with an interpulse interval T and varying response 
7. This form of F{t) compares release protocols with equal long-time-averaged probe 
release. For the parameters chosen, the constant release F{t) = 6{t) results in the cell 
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arriving at the near target. For pulsed release with interpulse interval T — 3, 7, the 
cell initially moves towards, and depending on the pulse intensity, may first reach the 
near target before eventually being pulled to the farther, stronger target. If the release 
frequency is increased even further, (interpulse interval T = 0.2, red dotted curve), 
the trajectories will again arrive at the closer, weaker target. Since the underlying 
processes are dissipative, at very high frequencies, the cell cannot respond fast enough 
to distinguish the pulses and a rapid succession of pulses is equivalent to an effective- 
amplitude, constant emission. 

In Fig. 4(c), we show a phase diagram indicating the regions of ^-T space in which 
we expect the cell to go to the near, weak target, or to the far, strong target. When 
the interpulse time T is small and the cell's response 7 is small, it will go to the far 
target. For 7 > 0.225, the near target is selected. These results are consistent with 
those shown in Fig. 4(a). When the interpulse interval T is large, the cell will reach 
the far target for any response strength 7. In this example, the cell will reach the far 
target for any value of T > 4. When T falls in the interval < T < 4, the cell will 
reach the far target when 7 is either small or large, but will reach the near target at 
intermediate values of 7. When the cell reaches the far target and 7 is large, initial 
chemoattractant pulses quickly pull the cell towards the near target. For several pulses, 
alternating chemoattractant waves from the far and near targets will pull the cell away 
from, then back towards, the near target. The cell appears to "bounce" around the near 
target. Eventually, a wave of chemoattractant from the far target will dislodge the cell 
from the near target, and the next wave of chemoattractant from the near target will be 
insufficient to return the cell to the near target. After this, each subsequent pulse brings 
the cell closer to the far target, until the cell reaches this target. The trajectory with 
parameters 7 = 1, T = 3 illustrates this qualitative behavior in Fig. 4(b). When the 
cell with small 7 reaches the far target, it does so without first reaching the near target. 
The trajectory with parameters 7 = 0.2, T = 7 in Fig. 4(b) illustrates this qualitative 
behavior. 

A more quantitative understanding of the target selection phenomenon can be found 
in the small 7 limit, where the cell does not move much under the influence of a single 
probe pulse (such as the trajectory corresponding to 7 = 0.2 and F{t) = 7^,- S{t — 7i) 
in Fig. 4(b)). If the cell does not move appreciably under the influence of a single probe 
pulse, we can use Eq. 7 to approximate the asymptotic distance traveled by the cell as 
a result of one single probe pulse released ai t — 0. If < 0, subsequent probe pulses 
will be released when the cell is closer to the left target, and the cell will eventually 
incrementally move leftward. If X^x, > 0, the cell will ultimately arrive at the right 
target. Therefore, X^ — defines an approximate boundary for selection between two 
targets in d = 1. 

Generalizing the phenomenon to d dimensions, we expect that for each release 
sequence F{t), there will be at least one d — 1 dimensional surface that separates 
trajectories that evolve to different targets starting from a given initial position on 
the d dimensional manifold. Figure 5 shows trajectories of cells searching and selecting 
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Figure 5. Target selection in c? = 2 depends on probe release protocol. The separatrix 
corresponding to 9{t) probe release (dotted line) divides the space such that trajectories 
(solid red lines) originating from points to the right of the dotted line arrive at 
the stronger target (big open circle) while those starting from the left arrive at the 
weaker target (small open circle). Amplitude-dependent target selection is shown by 
using F{t) = 0.1 X d{t) (black dashed line), leading to a different target than when 
F{t) — 9{t). Likewise, frequency-dependent target selection is demonstrated using 
F{t) — 10^j(5(i — lOi) (black solid line with a unit average probe release) where 
the cell selects the nearer target when started in a region that would normally lead 
to selection of the stronger, farther target when F{t) = 9{t). For comparison, the 
separatrix for continually constant probe release F computed by integrating along the 
ridge of the static field rib(r) is shown by the thin blue curve. Parameters used were 
D = ki — 1, k2 = 0.2, iia = fJ'b = 0.001 and 7 = 10. The numerical solution to Eqs. 
1-3 in d = 2 were found using the Peaceman-Rachford algorithm. 



between two targets in d = 2. The separatrix associated with the Heaviside probe 
release F{t) = 6{t) is indicated by the dotted curve that divides the space such that 
trajectories originating from the right/left of this hne are led to the stronger /weaker 
target respectively. For comparison, the separatrix for constant probe release for all 
times (corresponding to the 7 — )■ limit) computed by integrating along the ridge of 
the static field nf,(r) is shown by the thin blue curve, highlighting its sensitivity to F{t). 
We then start the cell within regions where one target is clearly selected when the probe 
chemical is released with F{t) = 0{t), but with either a diminished released rate or with 
pulsed release, such that the average rate of chemical released per unit time is 1. As in 
the d = 1 case we find that for Heaviside function release, a larger release rate favors 
the weaker, closer target, while a smaller release rate favors the stronger, farther target 
(thick dashed curve). In one dimension we find that slow pulsed release favors the far, 
strong target over the near, weak target. In contrast, for d = 2, we find that slow pulsed 
release favors the near target (solid black curve). This may be understood as follows. 
In d = 1, pulsed release causes the cell to quickly go towards the near target, but a 
subsequent wave of chemoattractant from the far target caused the cell to ultimately 
select the far target. In two dimensions, the radial divergence of the concentration fields 
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renders the signal from the far target insufficient to pull the cell from the weaker target. 
In higher dimensions, the radial spread is stronger, and the frequency and 7-dependent 
selection mechanism is significantly mitigated. 

While nontrivial target selection phenomena still arises in = 2, many cellular 
chemotactic responses, such as development, occur in d = 3 where the concentration 
fields spread significantly and attenuate signalling responses. However, in many 
situations, the medium in which the concentration fields Ua and rib diffuse is 
heterogeneous. A common biological example is the extracellular cellular matrix and the 
existence of intervening cells and tissues. Diffusive transport in a random medium can 
also be described by diffusion in a lower, effective dimension. Furthermore, media that is 
ramified can transport signals in a one-dimensional manner along the main percolating 
path once side branches have saturated with probe and/or chemoattractant. Therefore, 
signatures of our proposed mechanism may nonetheless arise in = 3. 

Although our results are based on a simplified model of response (Eq. 4) where 
U[nb{r,t)] oc nfe(r,t), and neglect noise of memory effects in 7, they provide clues 
into how more complex models might behave. For example, in the limit of a motility 
response that requires highly cooperative binding of multiple chemoattractants (large 
Hill coefficient a), the motion of the cell would be appreciable only when a very narrow 
range of gradients in n?, cross the cell. High cooperativity imposes a switching response 
which would lead to saltatory movements of the cell. While larger velocities can be 
attained (because VU[nb] can be larger is a is large) the time scale over which Vf/[nf,] 
is large is short. We expect the net effect to be quantitative and that most our qualitative 
results will persist. 

Delayed response, or directional memory, may arise from receptor adaptation and 
other intermediate steps in the cell signaling pathway [20] , and may give rise to behavior 
that is more difficult to intuit without numerical computations. However, simple 
conclusions in certain limits can be motivated by comparing the delay time, or directional 
correlation time with the arrival times of the gradients VC/fri;,] from the targets. For 
example, if directional memory is longer than the difference in times it takes for the Ub 
signal to propagate from the close and far targets, the cell would be less sensitive to the 
signal from the far target. This would have the effect of biasing selection of the near 
target, provided the release frequency is sufficiently high. 

One example of a biologically-motivated variant of our model worth further 
discussion is "autocrine" signaUing. Autocrine signalling is often relevant to bacterial 
aggregation, where the probe is identical to the chemoattractant {a — b). Within our 
mathematical framework, the autocrine mechanism is described by replacing nb{r,t) 
with na{r,t) in Eqs. 2 and 4, and setting fib = A*a Db = Da- This model is also 
novel in that its naive continuum limit does not map onto existing PDE descriptions 
such as Kcllcr-Segel type models. In the Keller-Segel model, the production of autocrine 
chemoattractant is simply proportional to the local cell density and is independent of 
the probing signal the cells receive. Our qualitative results for paracrine signalling also 
hold in the autocrine signalling problem, which can be understood by considering the 
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Ua profile generated by the cell and that generated by the target. In the autocrine case, 
when the cell responds to n^, its self-produced field is maximal near the point where 
it is released. Therefore, the self-produced tends to keep the cell from moving. When 
the gradient in generated from the target is finally felt, the cell will move towards 
the target; however, typically by this time, the self-generated Ua has diffused and/or 
dissipated away such that it's "restoring force" is weak. Therefore, we expect that all 
else being equal, the replacement of rib by ria, leads to an overall slightly suppressed 
chemotactic response. 

Finally, analysis of the effects of stochastic responses, implemented through additive 
Brownian terms to Eq. 3 or 4, or through a time random response coefficient 7(t), is 
beyond the scope of this study. However, it has been shown that in a two-dimensional 
autocrine mechanism, stochastic effects and chemical decay keep the cells diffusive [18], 
despite the tendency for all cells and targets to aggregate. Therefore, one might expect 
that target selection to be suppressed in the presence of sufficient noise. 

4. Summary & Conclusions 

In conclusion, we have proposed a model for dynamic, multistep chemotaxis that 
involves chemical communication between cell and targets. Our analysis shows that 
signalling agents can select among potential targets by controlling the amount of, and 
frequency at which probe chemical is released. Since our moving source problem is 
intrinsically nonlinear, we employed numerical calculations to provide evidence of a 
critical target-switching release amplitude, as well as a window of response strengths 7 
within which a cell chooses the nearer, weaker target over the farther, stronger target 
that is selected outside this window. This effect arises from a nonlinear interplay 
between diffusion, decay, and chemoattractant production and depends on the frequency 
of release. Moreover, we found that for the conditions studied, a constant probe 
release leads to acquisition of the weaker, nearer target, while low frequency spike 
release leads to acquisition of the stronger, farther target. At higher frequencies, 
the cell again approaches the weaker, nearer target since high frequency and constant 
release give rise to similar spatial-temporal probe distributions due to the averaging 
nature of the underlying diffusive physics. Our numerical experiments have shown that 
target selection is observed over a wide range of system parameters, suggesting the 
chemolocation mechanism may be common in Nature. 

Numerous variants of our basic model, such as mobile targets that sense 
probe chemical, stochastic effects, and delays in the signaling processes can be 
straightforwardly investigated. In the present work, we assumed the simplest 
deterministic, instantaneous response in order to avoid introducing excessive number 
of parameters. Nonetheless, we have argued that many extensions of our underlying 
model, except sufficiently large noise, are likely to qualitatively preserve the predicted 
results. 
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